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1 Introduction 

On April 2015, the J-PARG E56 (JSNS^: J-PARC Sterile Neutrino Search using neu¬ 
trinos from J-PARC Spallation Neutron Source) experiment officially obtained stage-1 
approval from J-PARC. We have since started to perform liquid scintillator R&D for 
improving energy resolution and fast neutron rejection. Also, we are studying Avalanche 
Photo-Diodes (SiPM) inside the liquid scintillator. In addition to the R&D work, a back¬ 
ground measurement for the proton beam bunch timing using a small liquid scintillator 
volume was planned, and the safety discussions for the measurement have been done. 
This report describes the status of the R&D work and the background measurements, in 
addition to the milestones required before stage-2 approval. 

2 Experimental milestones before stage-2 approval 

There are two main points to show the J-PARC PAC before stage-2 approval. 

(1) The concrete detector location should be considered including the physics potential 
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of the experiment. This was written in the minutes of the 19th PAC [I]. For this pur¬ 
pose, we have to discuss this issue in close coordination with the facility people, since the 
constraints from the facility are important to consider alongside the physics ones. 

(2) We assume a rejection factor of 100 for fast neutron events induced by cosmic rays. 
The technique using Cherenkov light at LSND already achieved this factor, however 
we have to show the performance of the liquid scintillator (and photo-detector) quan¬ 
titatively. Both the Cherenkov technique and the Pulse Shape Discrimination (PSD) 
technique are tested for this purpose. 

For the former, we would like to have a special joint discussion with PAC members 
and MLF facility people (the FIFC is a candidate). The discussion includes safety issues 
as well as facility maintenance. We will show the status of the R&D for the latter. 


3 R&D for the JSNS^ Detector 


3.1 Motivation 

As shown in the status report submitted to the PAC last November |2[ |3l 0], all 
backgrounds in the candidate detector location are manageable for performing the JSNS^ 
experiment. One crucial assumption is that the detector should have a rejection factor 
of 100 for fast neutron events induced by cosmic rays. Table shows a summary of the 
number of background and signal events in the JSNS^ (shown in reference |2]). Based 
on these assumptions, the background rate induced by cosmic rays is comparable to the 
other background rates. 



Contents 

events 

Signal 

sin‘^2e = 3.0 X IQ-^ 
Am^=2.5eR^ 

(Best fit values of MLF) 

480 

sin^2e = 3.0 X 10“^ 
Am^ = 1.2eV^ 

(Best fit values of LSND) 

342 

background 

z7g from iJi~ 

beam-associated fast n 
Cosmic-induced fast n 
Accidental coincidence 

237 

16 

<13 

37 

32 


Table 1: Number of events in 50 tons. IMW beam x 5 years is assumed. 


The factor of 100 fast neutron rejection capability was already achieved by the LSND 
experiment using Cherenkov light detection, since protons in the energy range of 20-60 
MeV (recoiled by fast neutrons induced by cosmic rays) cannot emit Cherenkov light. 
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However we have to show the performance of JSNS^ detector independently since the 
design is different from the LSND. 

We are also trying to establish a PSD technique, because it uses independent infor¬ 
mation from Cherenkov light: The output signals’ pulse shape from liquid scintillator 
are different between protons recoiled by fast neutrons and positrons from Inverse Beta 
Decay (IBD). If we can use both information from Cherenkov and PSD, the rejection 
power could be even stronger than what has been assumed. 

In the following subsections, we show the status of the R&D using small prototype 
devices. For the Cherenkov test, there is a ~10L prototype, and a O.IL level prototype 
is used for the PSD test. 

The possibility to use APD (SiPM) is also being considered, since the space of the 
detector is tightly restricted, thus SiPMs may be used for the veto region even though 
the photo-coverage is small. Therefore, its R&D status is also described briefly. 


3.2 PID Capability using Cherenkov Light in Liquid Scintillator 

The JSNS^ experiment plans to use linear alkylbenzene (LAB) as the base of the 
liquid scintillator. However, we have started to follow the measurements in reference [H] 
using Mineral oil -|- 0.03g/l b-PBD at first, since the LSND study will be an excellent 
reference for subsequent measurements. A cylindrical prototype detector with dimensions 
of 130mm (diameter) x 1000mm (height) was prepared for this test. Figure shows the 
measurement setup. The cylindrical prototype is filled with water or diluted scintillator, 
and a black sheet is lining inside the detector to avoid reflections of light. The motivation 
of this experiment is to measure the timing difference between Cherenkov light and scin¬ 
tillation light mainly. The Cherenkov light has a faster emission timing than that of the 
scintillation light in general. Note that water does not emit scintillation light, therefore 
it provides a good cross check of the timing information. 

If the photo-detector receives the scintillation and Cherenkov light at the same time, 
the timing distribution should be like Fig. The solid line shows the scintillation light 
while the dashed line corresponds to the Cherenkov light. Here we assume that the time 
distribution of the scintillation light emission is the sum of a exponential for the fast time 
component and a (l-l-t/r)"^ empirical shape of the slow light, as written in the LSND 
NIM paper IB0 We also assume a time resolution of 0.5ns for the photo-detector (only 
for Cherenkov). Note that these time constants depend on the liquid scintillator being 
used. 


3.2.1 Cross-check with filling Water 

For the water case, there is no scintillation light, and we therefore did not see any 
light in the setup shown in the bottom of Fig. [^b). On the other hand, we can see 
Cherenkov light with the setup shown in the top of Fig. [^b). 

^ The fast time constant is 1.65 ns, and the slow time constant is 22.58 ns. The mixture ratio is 57% 
and 43% for each component. 
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(a) Setting for Cherenkov + scintillation (b) Setting for scintillation light only 

Figure 1: Schematics of the test setup, (a) shows the setting for the Cherenkov plus 
scintillation detection and (b) shows the setting for the scintillation light only detection. 
The cylindrical prototype with a size of 130mm (diameter) x 1000mm (height) is filled 
with water or diluted scintillator (blue parts of the schematics). Scintillation counters, 
which are shown as red boxes, are used to tag cosmic ray muons. Coincidence signals 
by the scintillation counters are used as the reference timing for the Cherenkov and 
scintillation light inside the cylinder. The effective area for each scintillator is 5 x 10 
cm^. The average muon path distance from the 2 inch PMT was kept as 65 cm to make 
the amount of scintillation light the same in both settings. 


Figure]^ shows the typical waveform timing distributions in an oscilloscope. We will 
use the oscilloscope data for the analysis since the cosmic ray trigger rate is very low. 
The two plots show typical timing distributions of the coincidence from the scintillation 
counters (top) and Cherenkov light (bottom). Note that the coincidence signal from 
the scintillation counters was delayed due to electronics. Also, this relative timing was 
smeared by the PMTs jitter (~1 ns), electronics (within 1 ns), and hit position of the 
cosmic rays (<< ns). This measurement is important for understanding the effect of the 
jitters. 

If the pulse height exceeds below -2000 counts beyond the ground line (e.g.: ~20000 
is the ground line, and the threshold is ~18000 in the bottom plot of Fig.|^, it is regarded 
as signal. -2000 corresponds to a threshold of about 1/3 photo-electrons (p.e.), therefore 
the threshold is low enough to search for the signal hits. 

Figure shows the relative timing of water Cherenkov signals with respect to the 
coincidence. The timing of the fastest point of each signal which exceeds 2000 counts is 
plotted in this figure. Out of 100 triggers, about half of them have signals beyond the 
threshold, and this rate is reasonable according to the MC simulatiorj^ As shown in the 
plot, the total jitter of the water Cherenkov signal (or the relative timing resolution from 
systematic effects) is ~ 1.8 ns. 

^ MC calculation shows that the Cherenkov light gives 0.5^0.8 p.e. / trigger. The uncertainty is due 
to the momentum and angle of cosmic rays. They are not included in the MC yet because the timing 
does not depend on those. 
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Figure 2: Timing distributions of the scintillation light (solid) and Cherenkov light 
(dashed). The horizontal axis corresponds to the time, unit is ns. 

3.2.2 Diluted Scintillator 

Hereafter we show the measurements using a diluted scintillator, which contains min¬ 
eral oil plus 0.03g/l of b-PBD. This diluted scintillator was used in the LSND experiment, 
and it can be used as the reference of the other measurements or for discussion of detector 
performance. 

Setting for scintillation light only detection (Fig. [^b)) 

The scintillation light with the setting shown in Fig. [^b) is discussed hrst. With 
this setting, there should be only scintillation light detection. Figure shows the typical 
signals using the diluted scintillator. Events with a slow scintillator timing signal and 
multiple light signal are shown. 

To analyse events with multiple hits, one hit is dehned as follows. 

• Points less than -2000 counts are searched for. 

• If there are more than two points, it is judged whether they are in sequence in time 
or not. If they are in sequence, they are regarded as one hit. 

• The fastest point within one hit is recorded as the “timing of the hit” 

Figure shows the number of hits per one event and the hit timing distribution. 
Clearly, the slow components from the scintillator can be seen. 

Setting to detect scintillation plus Cherenkov light (Fig. [^a)) 

An analysis with the setting shown in Fig. [^a) was performed as well. The procedure 
for this analysis is exactly the same as what was described before. Figure shows the 
results on the number of hits and timing of the hits. As expected, both the number of 
hits (left) and the number of hits around the fast timing (right plot - same timing as the 
expected Cherenkov timing), increased. 
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Figure 3: Typical timing distributions of the coincidence from the scintillation counters 
(top) and Cherenkov light (bottom) measured by the oscilloscope. The horizontal axis 
corresponds to the time. One count is 2 ns. The vertical axis of the plots is the pulse 
height counts. 1000 counts corresponds to ImV in the bottom plot for this case. 



Figure 4: Timing of the Cherenkov light. The horizontal axis is time counts, and 1 
count corresponds to 2 ns. With respect to the coincidence signal, the total jitter is 
<0.9 counts (1.8 ns) from the ht result as shown in the bottom of the fit parameters. 
(0.8998±0.1302) 
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(a) Typical event 1 (slow timing) (b) Typical event 2 (multiple light) 

Figure 5: Typical signals using diluted scintillator. The horizontal axis is time counts, 
and 1 count corresponds to 2 ns. (a) We see a slow timing hit, and (b) multiple hits are 
seen. Note that the Cherenkov signal is seen around ~24982 counts on the horizontal 
axis as shown in Fig. 



(a) Nhit distribution (b) Hits timing distribution 

Figure 6: (a) Number of hits in one trigger, (b) Hit timing distribution. The horizontal 

axis is time counts, and 1 count corresponds to 2 ns. Comparing with Fig. a clear tail 
(slow component light) from diluted scintillator is present. 


















(a) Nhit distribution (b) Timing distribution 

Figure 7: (a) Number of hits in one trigger, (b) Timing distribution of the hits (for 

Cherenkov+scintillation). The horizontal axis is time counts, and 1 count corresponds to 
2 ns. 

Comparison between scintillation only and Cherenkov+scintillation light 

Figure]^ shows the timing comparison between the scintillation light only case and the 
Cherenkov+scintillation light case. The horizontal axis is changed to ns to understand 
it better (starting point is still arbitrary). Also both data sets are normalized using 
the number of triggers. Two things are easily understood: (A) Cherenkov timing is 
faster than scintillation light as expected. (B) The amount of the Cherenkov light is 
comparable to the scintillation light in this scintillator condition. For the future work for 
item (A), we have to compare this result with Fig. [^quantitatively (with smearing 1.8 ns 
at least), since one may reduce total jitter. For the future work of item (B), to compare 
the amount of the light from Cherenkov and scintillation, a detailed MC simulation is 
necessary, especially for the Cherenkov part since cosmic rays have momentum and zenith 
angle distributions. 


3.2.3 PID Summary 

We checked the performance of a LSND type liquid scintillator directly. The Cherenkov 
light is useful information in this diluted scintillator condition as expected. However, this 
statement is still qualitative, and we need a more quantitative statement in combination 
with the MC simulation. 

Three future directions should be examined with this R&D. (1) variation of the 
concentration of the scintillator components for diluted scintillator, (2) similar study 
using linear alkylbenzene (LAB), and (3) test beam with protons. All will be done in 
order with priorities. 

With these inputs, it is important to estimate the realistic PID capability using MC 
simulation with a real size detector. This will be done within about 1.5 years. 
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Figure 8: Timing comparison between scintillator only (cross) and Cherenkov + scin¬ 
tillator (histogram). Cherenkov light is concentrated in the fast timing, and the yield of 
Cherenkov is similar to the scintillation light. The horizontal axis is changed to ns in 
order to improve understanding (starting point is still arbitrary). 

3.3 Pulse Shape Discrimination (PSD) of Liquid Scintillator 


3.3.1 R&D of Gd loaded liquid scintillator(GdLS) 

The Japanese group tried to load Gd to LAB based liquid scintillator following the 
Daya-Bay experiment’s paper [7j. The JSNS^ experiment will use this liquid type and 
the test with a small size detector will be done in Japan. In order to load Gd in LAB, it is 
necessary to use a Gd-carboxylate complex formed by a mixture of gadolinium chloride, 
neutralization solution with TMHA, and ammonium hydroxide. A sample of the Gd- 
carboxylate complex was made and used to produce the Daya Bay type liquid scintillator 
with 0.1 w% of Gd concentration. 

In order to check whether Gd was successfully loaded in the LAB, the thermal neutron 
capture time was measured since there is a correlation between the neutron capture time 
and Gd concentration. 

The measurement setup is shown in Fig. 

The volume of the sample is lOL contained in an acrylic box. The scintillation light is 
viewed by two 8 inches PMTs and the signals are taken with a flash ADG module (CAEN 
V1721, 500MS/s, 8bits). An ^^^Am®Be source, which emits both a 4.4 MeV gamma and 
neutron below 11 MeV, is set at the center of the acrylic box. The gamma and protons 
recoiled by the neutron are detected as the prompt signal, and gammas from thermal 
neutron capture on Gd are detected as delayed signal after about 30/isec of mean time. 


Figure 10 shows the time difference (At) distribution between the prompt and delayed 
signals, which includes samples of accidental coincidences. Therefore, it was fit with 
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Setup 



Figure 9: Picture of the setup for measuring the thermal neutron capture time with the 
Daya Bay type GdLS. 


a sum of an exponential and flat function. Then the thermal neutron capture time is 
estimated as 29.8±0.2 /isec, which can be compared with the expected time of 28 fisec. 
Considering the lit result, it was confirmed that about 0.09 w% of the Gd concentration 
was loaded in LAB. 



Figure 10: Distribution of thermal neutron capture time using the Daya Bay type GdLS. 
The blue line shows the flt result with a sum of exponential and flat functions. 


3.3.2 R&D of GdLS with capability of n/y pulse shape discrimination (PSD) 


Measurement of the PSD capability for Daya Bay type of GdLS 

We measured the PSD capability for Daya Bay type scintillator(DBLS) without Gd 
(only LAB, PPG and bisMSB). The neutrino energy in JSNS^ is a few tens MeV. The 
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PSD capability increases with the neutrino energy, so if the DELS shows enough cosmic 
induced fast neutrons rejection power while maintaining high signal efficiency in the 
energy range around a few tens MeV, we can use DELS without any PSD capability 
improvement. 

Figureshows the setup of the PSD capability measurement. The target scintillator 


Setup for PSD measurement 


LSforTOF 
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Figure 11: Scheme of the PSD capability measurement. 


is contained in a 100 mL vial and another liquid scintillator volume, with 1 L, where the 
^^^Cf source is attached. The ^^^Cf source emits both 7 and neutron with an energy of 
few MeV, and the coincidence signals of the scintillators are detected for the PID. The 
distance between the target scintillator and the ^^^Cf source is 40cm. Using the hit time 
difference between both scintillators, the neutron events can be distinguished from the 
gamma events, because the time of flight(TOF) of the neutron is different from that of 
gammas, as shown in the vertical axis of Fig. [T^ 

Waveforms of both scintillators are taken with a flash ADC module (CAEN V1730, 
500MS/s, Mbits). To evaluate the PSD capability, the ratio of the tail integrated charge 
of the waveform to the total integrated charge is used (TailQ/TotalQ variable). In this 
analysis, events around 500 photoelectrons are used for the evaluation of PSD capability 
since this number of photoelectrons corresponds to the mean energy of neutrons / 7 s of 
the ^^^Cf source with consideration of the light acceptance of the 2 inches PMT attached 
to the vial vessel. On the other hand, this number of photoelectrons corresponds to 3.5 
MeV in the JSNS^ experiment considering the acceptance difference of the detectors. 
Neutrino energy is a few tens MeV in the JSNS^ experiment, so this analysis gives a 
very conservative evaluation. To evaluate the PSD capability in the actual energy range, 
higher energy neutron and electron sources are necessary. 

Figure [I^ shows the correlation between the TailQ/TotalQ variable and the detected 
number of photoelectrons for the DELS. In the plot, there are two bands along the 
horizontal axis, which show the separation between gammas and neutrons. The left 
plot of Fig. 13 shows the correlation between the TOF and the TailQ/TotalQ variable 
around 500 photoelectrons in Fig. and the right plot shows the ID histograms of 
the TailQ/TotalQ distributions for gammas and neutrons. Using the TOF cut (orange 
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Figure 12: Correlation between the TailQ/TotalQ variable and the detected number of 
photoelectrons for the DELS. 



Figure 13: A correlation between the TOF and the TailQ/TotalQ variable(left plot) 
around 500 photoelectrons in Fig. and ID histograms of the TailQ/TotalQ distribu¬ 
tions for gammas and neutrons (right plot) 
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line in the left plot), the neutron events can be distinguished from the gamma events. 
Each TailQ/TotalQ distribution after applying the TOF cut was fit with a Gaussian 
function, and the rejection power of neutrons and the detection efficiency of neutrino 
events (gamma) can be estimated by using the fit results. Figure [l4| shows the correlation 
between the detection efficiency of gammas and the neutron mis-ID fraction. At a factor of 



Figure 14: Correlation between the detection efficiency of gammas and neutron Mis-ID 
fraction. 


100 rejection power for neutrons (neutron Mis-ID fraction ~ 0.01), the detection efficiency 
of gammas for the DBFS is 0.85. The high detection efficiency using the DBFS is shown 
around 3.5 MeV equivalent already, so an efficiency increase is expected in the energy 
range of a few tens MeV due to photo-statistics. 


Future plan 

It is important to measure the PSD capability in the actual JSNS^ energy range of 
a few tens MeV. Neutrons and electrons created by an accelerator can be considered for 
this test. 

In addition to testing with a small size detector, it is essential to estimate the PSD 
capability of the real size detector using MC simulations. Because of the distortion of 
the hit time distribution (or sum of waveforms) due to scintillation light attenuation, the 
light reflection by the tank and mis-reconstruction of the vertex, a large detector has 
worse PSD capability. 

We need to design a detector which can precisely measure the time distribution of 
scintillation fight emission event by event for high PSD capability. We are considering 
the following crucial items; 

(A) Develop a tank with less light reflection 

(B) Flash ADC module with high sampling rate(>500MS/s) and number of bits(>10bits). 

(C) Optimization of number of photoelectrons (light yield of the GdFS). 

(D) Develop a good vertex reconstruction algorithm. 
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3.4 APD (SiPM) 

The current design of the JSNS^ detector is based on technologies that have suc¬ 
ceeded in other experiments. The geometric size of the JSNS^ detector, however, will 
be comparatively rednced to adapt to the MLF facility. With respect to the optical de¬ 
tection, large photomnltipliers (10 inch) will be nsed as photon-detectors. In this case, 
large-sized PMT could affect the design of a snb-detector, snch as a Veto connter. This 
inflnence becomes large in the compact detector system, snch as with JSNS^. Therefore, 
we have started to stndy the possibility of using APD (avalanche photodiode) technology 
as a snbstitnte of the photomultiplier. It is one possibility as a supplement to the JSNS^ 
detector design. 

For the JSNS^ test, Hamamatsu photonics provided some engineering samples of 
MPPC (multi-pixel photon counter), APD production. We have started the verification 
of the MPPC in liquid scintillator for long-term usage. At present, three types of MPPC 
were tested (listed in Table. [^. 


Table 2 : API 

D list 

Model 

Package Type 

Pixel 

Window material 

S-10362-11-025U 

Metal 

1600 

Borosilicate glass 

S-10362-33-050U 

Ceramic 

400 

Silicon 

S-10362-11-100U 

SMD 

100 

Epoxy resin 


In the JSNS^ experiment, liquid scintillator(LAB) is chemically active. A MPPC 
may therefore be damaged by the liquid scintillator. We are studying potential damage 
to the window material due to the package type. 

and continue to observe 


For the first test, we used the mineral oil in section 3.2 


an influence on the window of MPPC in liquid scintillator. The silicon window of the 
ceramic package was deformed and become cloudy in a day. 

In the future, we will evaluate the influence with a term of a month or year. We 
will check the deformation of the window and gain of the MPPCs. We have a plan for 
development of a MPPC in the liquid scintillator in collaboration with the manufacturer. 


4 Background Measurement on the Proton Bunch 
Timing using 1.6L detector 

As shown in the previous PAC [2], there is a significant neutral particle flux that is 
observed in the proton bunch timing window, while this background is manageable outside 
of this window for the JSNS^. An additional particle ID measurement (neutrons or qs) 
for the particles on the proton bunch timing using a small size (1.6L) liquid scintillator 
detector was planned at the MLF 3rd floor. We assnmed that all particles are nentrons 
in the previons status report, but, if there is some 7 fraction, the nnmber of backgronnds 
for the delayed signal will be smaller. 

Unfortnnately, there have been two large incidents (fire from muon site and mer¬ 
cury target cooling issue) at the MLF. Therefore, the measurement was not done yet, 
but the safety discussion was successfully done carefully between facility people and the 
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experiment. Hopefully, this measurement will be done starting this autumn. 


5 Summary and Plan 

In this report, we show the status of our R&D work, which is going well. How¬ 
ever, the concrete detector performance with quantitative numbers with the realistic MC 
simulation should be shown before the stage-2 approval. 

In addition to this R&D work, a realistic detector location will be investigated in 
consideration of the the facility constraints. 
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